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Abstract 

We study proton decay in a supersymmetric SO(10) gauge theory in six dimensions 
compactified on an orbifold. The dimension-5 proton decay operators are forbid- 
den by R-symmetry, whereas the dimension-6 operators are enhanced due to the 
presence of KK towers. Three sequential quark-lepton families are localised at the 
three orbifold fixed points, where SO(10) is broken to its three GUT subgroups. 
The physical quarks and leptons are mixtures of these brane states and additional 
bulk zero modes. This leads to a characteristic pattern of branching ratios in proton 
decay, in particular the suppression of the p — > fi + K° mode. 



1 Introduction 



The most striking consequence of Grand Unified Theories is proton decay [1-6]. It has 
been predicted and sought for more than 30 years, and its absence constitutes a very 
strong constraint on any realistic GUT model. 

Recently, the discussion of proton decay has been revitalised on two different fronts. 
On the experimental side, the bounds coming from the SuperKamiokande experiment 
have reached r{p — > e + 7r°) > 5.3 x 10 33 yrs [7] and r(p —> K + u) > 1.9 x 10 33 yrs [8]. On 
the theoretical side, these bounds have motivated new detailed studies of proton de- 
cay via dimension-5 operators in supersymmetric SU(5) models [9-13]. These analyses 
showed that the present bounds disfavour this class of models, although the theoretical 
predictions strongly depend on the flavour structure assumed [12,13]. 

Dimension-6 operators are less dangerous and have not drawn so much attention 
in the literature. It is, however, well known that such operators depend on the flavour 
structure as well [14-16], and it has been realised that the observed leptonic mixing can 
have a strong effect on the proton decay branching ratios [17]. Recently, dimension-6 
operators have been studied in flipped SU(5) [18] and SO(10) [17,19] GUTs. 

The interest in dimension-6 operators has been renewed by GUT models in higher 
dimensions, with symmetry breaking via orbifolds. Here dimension-5 proton decay is 
naturally absent [20,21]. In such models though, the dimension-6 operators are enhanced 
compared to the usual 4D case, due to the lower mass scale for the heavy particles 
mediating the decay and the presence of Kaluza-Klein towers of such states. Furthermore, 
the proton decay branching ratios depend on the localization of quarks and leptons [22]. 

The goal of the present paper is to study proton decay via dimension-6 operators 
in a specific SO(10) orbifold GUT model in 6D, where the complete Standard Model 
flavour structure can be reproduced via mixing of brane states with bulk split multiplets 
[23]. This kind of model can naturally arise in compactifications of the heterotic string 
[24,25]. Quark and lepton mass matrices are approximately of lopsided type, with some 
characteristic modifications compared to the 4D lopsided picture. The up quark of the 
first generation is located on a brane where the bulk gauge group SO(10) is broken to 
SU(5)xU(l), and therefore SU(5) gauge boson exchange gives the main contribution to 
proton decay. For this model, we will calculate the total rate and the branching ratios of 
proton decay and compare the results with those in a 4D SU(5) models with U(l) family 
symmetry. 

The paper is organised as follows: in Section we discuss the flavour structure 
of the 6D orbifold GUT model and evaluate the mixing matrices needed to calculate 
proton decay widths. Section El deals with proton decay via dimension-6 operators in the 
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Figure 1: The three SO(10) subgroups at the corresponding fixed points (branes) of the orbifold 
T 2 /Z 2 x Z' 2 x Z'i. 

usual 4D well as the 6D model for which the sum over the Kaluza-Klein tower 

is performed. Finally, in Section |U we discuss the effects of the flavour structure on the 
proton decay branching ratios in the different models. Conclusions are given in Section^ 



We start from an SO(10) gauge theory in 6D with N = 1 supersymmetry compactified on 
the orbifold T 2 /(Z£xZ% s xZ% G ) [26,27]. The theory has four fixed points, 0„ PS , GG and 
Ofl, located at the four corners of a 'pillow' corresponding to the two compact dimensions 
(cf. Fig.HJ. At O u only SUSY is broken, whereas at the other fixed points, PS , GG and 
O fl , also SO(10) is broken to its three GUT subgroups G PS = SU(4) x SU(2) x SU(2) [1], 
G GG = SU(5) x U(l) x [2] and flipped SU(5) [5,6], G fl = SU(5)' x U(l)', respectively. The 
intersection of all these GUT groups yields the standard model group with an additional 
U(l) factor, G SM / = SU(3) x SU(2) x U(l) y x U(l) y/ , as unbroken gauge symmetry below 
the compactification scale. 

The field content of the theory is strongly constrained by imposing the cancellation of 
irreducible bulk and brane anomalies [28,29]. We study the model proposed in [23], con- 
taining three 16-plets ijji, i — 1, . . . , 3, as brane fields, and six 10-plets, Hi, . . . , H 6 , and 
four 16-plets, $, $ c , (ft, 4> c , as bulk hypermultiplets. Vacuum expectation values of $ and 
$ c break the surviving U(1) B _ L . The electroweak gauge group is broken by expectation 
values of the anti-doublet and doublet contained in Hi and H 2 . 

We choose the parities of H 5 , H 6 and <f), cj) c such that their zero modes are 



These zero modes act as a (vectorial) fourth generation of down quarks and leptons and 
mix with the three generations of brane fields. We allocate the three sequential 16-plets 
to the three branes where SO(10) is broken to its three GUT subgroups, and place ipi at 
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GG , ^2 at Ofl and ^3 at PS . The three 'families' are then separated by distances large 
compared to the cutoff scale M*. Hence, they can only have diagonal Yukawa couplings 
with the bulk Higgs fields. Direct mixings are exponentially suppressed. The brane fields, 
however, can mix with the bulk zero modes for which we expect no suppression. These 
mixings take place only among left-handed leptons and right-handed down-quarks. This 
leads to a characteristic pattern of mass matrices. 

As described in [23], after B — L breaking at the scale and electroweak symmetry 
breaking via V\ = (Hf), v 2 = (H2), the mass terms assume the characteristic form 



W = d a m a al3 d c p + e c a m e al3 ep + n c a m®pVp + i*im"< 



c N r 



(2) 



4, include also the bulk states with 714 



^4, while 



where the Greek indices, a = 1 
Latin indices only run over brane fields, % = 1 . . . 3. Note that there is no Majorana mass 
for n\ = v\ which originates from the 6D hypermultiplet. m u and m N are diagonal 3x3 
matrices (tan/3 = i^/fi). 
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whereas m d , m e and m D are 4x4 matrices with the common structure, 
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Here fii,Jli = 0{v\) and M, = O(gut)- The diagonal and the off-diagonal elements of 
these matrices satisfy several relations due to the underlying GUT symmetries [23]. The 
hypothesis of a universal strength of Yukawa couplings at each fixpoint leads to the iden- 
tification of the diagonal and off-diagonal elements of m u / tan/3, m d , m e and m D / tan/3 
up to coefficients of order one. This implies an approximate top-bottom unification with 
large tan P and a parametrization of quark and lepton mass hierarchies in terms of the 
six parameters yUi,/i 2 ,yU 3 and /Ii,^,^. 

The crucial feature of the matrices m d , m e and m D are the mixings between the six 
brane states and the two bulk states. The first three rows of the matrices are proportional 
to the electroweak scale. The corresponding Yukawa couplings have to be hierarchical in 
order to obtain a realistic spectrum of quark and lepton masses. In a complete theory, 
this hierarchy may be due to the different location of the fixpoints on the orbifold. The 
fourth row, M a . is of order the unification scale and, as we assume, non-hierarchical. 



The mass matrix m (cf. Eq. (pj}) can be diagonalised by the unitary transformation 

m = U 4 U 3 DVM , (5) 

where the matrices U4, V4 single out the heavy mass eigenstate, that can then be inte- 
grated out, while the matrices U3, V3 act only on the SM flavour indices and perform the 
final diagonalisation in the remaining 3x3 subspace. t/4 and V4 are given by (neglecting 
phases) 
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M% + Mj. In general, V4 contains large mixings, 



where M = and M aP = 

while U4 is approximately the unit matrix up to corrections 0(v/M). U3 and V3 are the 
matrices that diagonalise 
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where 



m 



^2(^4)2^ + AA2<y 4 )4j 
\ /^3(V4)3j +A*3(Vk)4j / 



(9) 



Clearly, they have only a non-trivial 3x3 part, 

U 3 = 





CKM " 







Va 




(10) 



Notice that the rows of m scale each like /ij, /i;. Hence, for hierarchical parameters < 
^2,jj<2 < A*3)A*3 we obtain a structure familiar from lopsided fermion mass models. 
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As discussed in [23], we can choose the parameters in such a way to give a consistent 
quark mass pattern and CKM matrix, in particular 



Ail : A*2 : A*3 ~ ™u '■ m c : m t 



(11) 



and (cf. Eq. (jTTjl ) 

p 3 ~ m 6 , ju 2 : yU3 ~ m s : m b . (12) 

The CKM matrix in t/3 arises in a natural way as well. Setting the remaining parameter 
/xi to give 



/'1 



/'2 



0.2 



the other matrix elements are determined by the quark masses [23], 



H2 m s 

A*3 Wfe 



~ ^ ~ — ~ 2 x 10' 



Kb ~^ = e c ^^4x 10- 3 



(13) 



(14) 



/i 3 m 6 

Within the accuracy of our approach this is consistent with the analysis from weak 
decays [30]. 

The matrix V takes a relatively simple form in the case when /i lj2 = and the mass 
matrix in Eq. (JHJ) has one zero eigenvalue. We then obtain 
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up to a two-dimensional mixing matrix for the second and third generation, parameter- 
ized by a small angle Or, 



e, 



f4 + f4 „ 1 

Ms 



(16) 



where we have defined 



Ml \ 9 / Ml \ „ M 3 M A 



(17) 



The case of small /ii j2 limit is actually of physical relevance, since it gives for the 
down quark 



m d n 2 Hi n m c m b 

~ — — ~ H r — 0.03 
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consistent with data [31]. Since ^2/1^2 ~ fn c mb/(m t m s ) ~ 0.1, the corrections to the 
matrix (fT5|) are small. 

The charged lepton mass matrix m e has the same structure as the down-quark mass 
matrix, but there the large mixings are between the 'left-handed' states ej. Experimen- 
tal data require for the largest eigenvalue m T ~ m b , whereas the second and the third 
eigenvalue have to satisfy the relations m M ~ 3m s and m e ~ 0.2m^, respectively. This 
is consistent with our identification of m d and m e up to coefficients of order one un- 
less the relevant parameters are fixed by GUT relations. The comparision between the 
expression (fT%l) for and the corresponding one for m e suggests locating the second 
family on the flipped SU(5) brane. The successful relation for the light neutrino mass, 
m 3 ~ m^/M 3 ~ m^M^/vjf ~ 0.01 eV, requires the third family to be located on the PS 
brane. With the first family on the GG brane, m d and m e are determined by the param- 
eters \x\ and /if, which are not related by a flipped SU(5) mass relation. 

3 Proton decay via dimension-6 operators 
3.1 Effective SU(5) operators in 4D models 

Dimension-6 proton decay in the SU(5) model is mediated by the exchange of the X 
and Y leptoquark gauge bosons [16]. Contrary to the dimension-5 operator, it does not 
involve any dressing through supersymmetric partners and therefore it is not sensitive to 
the supersymmetry breaking scale (except for the weak dependence of the GUT scale on 
the superparticle mass spectrum). The effective vertex is obtained by simply integrating 
out the heavy gauge bosons. 

The couplings of the SU(5) representations 5* and 10 with the SU(5) gauge bosons 
are given by their kinetic terms, 



where g§ is the SU(5) gauge coupling; the Latin indices count the generations, z = 1 ... 3 
those of 10, k = 1 ... 3 those of 5*. 

We now express the SU(5) representations in terms of SM fields with X = (X, Y) 
being a (3*, 2, 5/6) representation of SU(3) x SU(2) x U(l). This yields the baryon and 




(19) 



which include 



C = z'-y= Al [2 tr (lOi 7 At T a 10i) + 5£ 7 ^(T a ) T 5£] + h.c. , 



(20) 
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lepton number violating terms 

C = -i^Xa» [e a/37 Qp ti Y u° ti + ~ c i Y Q a ,i - d? aM 7" L k ] + h.c. , (21) 

where Greek indices denote the colour degrees of freedom and the SU(2) indices have 
been suppressed. Note that the first two terms come from the 10 representation, the last 
one from the 5*. 

Integrating out the heavy gauge bosons with masses M%, we get the effective oper- 
ators relevant for proton decay 

g 2 — — — 

£ett = ~1^2 ^ QPJ [ e J Ta« Ql,j ~ d %k Jfi L k ] + h.c. . (22) 

With Fierz reordering, one can write the operators as 

£eff = ~Tf2 6 <x01 [ eC j uC a,i Qp,i Qj,j ~ d C a ,kU C /3,i Q-y,i L k ] + h.C. . (23) 

In the case of flipped SU(5) the X' = (X', Y') bosons form a (3*, 2, —1/6) represen- 
tation of SU(3) x SU(2) x U(l) with couplings 

C = -i-^K, Qf3 d c y - lF a Y L] + h.c. . (24) 

Contrary to SU(5), there is only a single baryon and lepton number violating operator 
in flipped SU(5), 

g 2 — — 

^eff = J7T2 e a/3^d c a!k U c p ti Q 7 ,i L k + h.C. . (25) 
M x 

3.2 Effective operators in 6D models 

In the orbifold model described above, the up-type quarks are localized at one fixed point 
each, in particular the up quark is located at the Georgi-Glashow one. It is therefore clear 
that dimension-6 proton decay can arise via the exchange of the SU(5) X and Y bosons 
as in the traditional 4D picture. There are though two important differences in the 6D 
case, as we will see in the following. 

First we have to take into account the presence of not only one X gauge boson, but 
of a Kaluza-Klein (KK) tower with masses given by 

(2n+l) 2 (2m) 2 

?2 P2~ 

for n, m — 0, 1, 2, . . .. The lowest possible mass is M#(0, 0) = I/R5, as given by the SU(5) 
breaking parity. Note that if we define the 4D gauge coupling as the effective coupling of 
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M 2 X (n, m) = V ; + y —f- , (26) 



the zero modes, the KK modes interact more strongly by a factor a/2 due to their bulk 
normalization. 

To obtain the low energy effective operator, we then have to sum over the KK modes, 

1 00 1 

2 £ 



2 E T • ( 27 ) 

n ^ (2n+l) 2 + §(2m) 2 



Taking formally the limit Rq/R^ — > 0, we regain the 5D result [22], 



(Mf) 2 ~(2n+l) 



The double sum in Eq. lj27Jl is logarithmically divergent. Since our model is valid only 
below the scale M*, where it becomes strongly coupled and also 6D gravity corrections 
are no longer negligible, we restrict the sum to masses M x (n,m) < M*. One easily finds 

K R 5 R 6 fin (M*R 5 ) + C (^f] + O ) . (29) 



(Mf) 2 4 °V v J \RdJ V^5/ 6 M, 

Note that in the logarithm the smallest KK mass appears, -R5 = 1/Mx(0, 0). The depen- 
dence of 1/(M|F) 2 on the cutoff M* has to disappear once the model is embedded in a 
more fundamental theory. In the symmetric case, R 5 = R 6 = 1/M C , one finds C(l) ~ 2.3 
and the expression (|29| simplifies to 

- UMA + 2. 3 ). (30) 



(Mf) 2 4M 2 V \M C 

Numerically, this agrees with the explicit sum over the KK masses within 1% for 
M*/M c = 10... 50, which is the relevant range for the ratio of cutoff and compacti- 
fication scales in 6D. 

The second, most important difference of 6D models compared to 4D models is the 
non-universal coupling of the X gauge bosons. In fact, due to the parities and the SO(10) 
breaking pattern, their wavefunctions must vanish on the fixed points with broken SU(5) 
symmetry, PS and Ofj, and therefore no coupling arises via the kinetic term with the 
charm and top quark or to the brane states d\ 3,^2,3- We also have couplings to the bulk 
states g?4, g?4 and However, due to the embedding of the zero modes in full SU(5) 

multiplets together with massive KK modes, i.e. (g^,/^), (g? 4 , L4), (DIJ4) and (D4JI), 
the charged current interaction always mixes the light states with the heavy ones, and 
it is therefore irrelevant for the low energy process of proton decay [22]. So the kinetic 
coupling in Eq. (|23|) arises only for a single flavour eigenstate, not for all flavours as in 
the usual 4D case. 
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3.3 Corrections from derivative brane operators 



Apart from the kinetic term couplings, at any brane additional couplings can arise con- 
taining derivatives along the extra dimensions of the locally vanishing gauge bosons. Such 
operators are a 6D generalization of the 5D derivative operators discussed in [22,32], 

C d = V 5 t (z) ! $1 {V 5/6 e 2V ) $ 2 + h.c. . (31) 

fixed point, M * J ^ 2 

Here c^ 6 are unknown brane coefficients, V 5 / e = d^/Q + iA 5 /Q are the covariant derivatives 
in the extra dimensions and $j are any two different, locally non- vanishing fields in group 
representations which form a singlet together with generators of the broken symmetries. 
For unbroken symmetries the chiral superfields A 5 / 6 vanish at the fixpoints. 

These supersymmetric terms produce on the GG brane couplings with the flipped 
SU(5) leptoquark gauge bosons X', whose derivatives do not vanish on that brane. On 
the flipped SU(5) brane, there are couplings containing the derivative of the X gauge 
bosons, and on the Pati-Salam brane there are derivative couplings with both X' and X. 
Due to these additional vertices, three different classes of operators can arise: 

• operators coming from X' exchange on the GG brane: these involve two derivative 
vertices and can produce additional contributions to the effective operator 

W k ¥ l Q 1 L j (32) 

with k,j = 1,4; we will discuss their contribution below; 

• operators coming from A'-exchange on the flipped SU(5) brane or X,X' exchange 
on the Pati-Salam brane: these usually involve the charm and top quark instead of 
the up quark, and they are therefore irrelevant for proton decay; 

• interbrane operators from both the exchanges of X and X' gauge bosons: they can 
involve either one or two derivative vertices and generate mixed flavour operators 
of the type 

t C j u\ U\ d k — d c k u\ Q\Lj k, j = 2,3,4 X exchange, (33) 

df«$ d 2 u k k = 1,4 X' exchange GG— fl, (34) 

u? di v k j, I = 3, 4 ; k = 1, 4 X' exchange GG— PS. (35) 

Apart from the last term, they are usually suppressed compared to the kinetic term 
operators by a factor M c /M* due to the different parities of the vertices. 
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To estimate the effect of these additional operators, which introduce a dependence 
on the 0(1) coefficients c 5 / 6 , let us now consider the KK summations with one or two 
derivative vertices. We restrict ourselves here to the case R 5 = R 6 = 1/M C and univer- 
sal coefficients c 5 / 6 at the different fixpoints. Note that even if suppressed by M*, these 
operators can be as important as the usual ones, since the derivative enhances the di- 
vergence of the KK summation, which compensates the suppression. For example, from 
the exchange of the X' bosons on the GG brane, we obtain the sum 

1 2 ^ |c 5 (2n + l) + c 6 (2m + l)| 2 

(Mf)L. (2n + l) 2 + (2m + l) 2 ' W 

which is quadratically divergent. Using again the cutoff Mx(n,m) < M*, we obtain 

* ' " ' + (37) 









- N 



{Mf)\ 16M C 2 V 7T 10 bi \M, 



to be compared with the 'canonical' term Eq. (|HUj) . A similar result is obtained for the 
exchange of X' gauge bosons between the GG and the PS brane, affecting only the decay 
into neutrinos. 

The exchange of gauge bosons between different branes involving only one single 
derivative operator are less dangerous since the propagator gives a factor (— l) n account- 
ing for the different parities of the derivative and the kinetic term vertices. Therefore 
the KK summations are in general suppressed by a factor M c jM*. The exchange of X 
bosons between the GG and the flipped SU(5) brane gives 

1 _ 2 ^ n c 5 (2n + l) + c 6 (2m) 



( M x)l,o. M cM*^ K ' (2n + l) 2 + (2m) 2 ' 

This sum is only logarithmically divergent thanks to the alternating signs. Using again 
the cutoff Mx(n,m) < M*, one finds 

1 W£Wl)Y (39) 



{Mf )l 2.\I,.\I, V \M t 

which is suppressed compared to the 'canonical' term by a factor M c /M*. The same 
result is obtained for the X' exchange between the GG and the flipped SU(5) branes. 

Note that in principle operators with a higher number of derivatives can also be 
present, which contribute at the same level as the single derivative ones because the 
divergence of the KK summation compensates the suppression by powers of M* . We will 
assume that their contribution is small. 

Regarding the N=2 scalar superpartners of the SU(5) gauge bosons, X^, the brane 
terms in Eq. (j3*T|) give rise to couplings with fermion kinetic terms which can only produce 
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corrections of order (mp/M*) 2 . The derivatives of the X 5j6 bosons do not couple to fermion 
pairs and are therefore irrelevant. 

Finally, we emphasize that the position of the lightest quark generation is crucial 
for the discussion of proton decay. For instance, if the up quark were located on the 
Pati-Salam brane, the dimension-6 operator coming from the kinetic terms would be 
absent since both the X bosons of SU(5) and the X' bosons of flipped SU(5) vanish 
there. In principle, this gives us a means to avoid the 'canonical' dimension-6 operators 
completely, leaving the derivative couplings as dominant contributions. 



4 Flavour structure and branching ratios 

4.1 Flavour mixing in 6D versus 4D GUT models 

Proton decay involves only the light quark states and the operators containing the com- 
binations uud and udd. Therefore we have to rotate the weak eigenstates into the mass 
eigenstates and single out the contributions for the lightest generation. Without loss of 
generality and for future convenience, we can start in the basis where the up-quark mass 
matrix is diagonal. Then the down quark and lepton mass matrices are not diagonal, in 
general, but can be diagonalised by unitary transformations, 

d L = U d L d' L , e L = U e L e' L , v L = U v L v' L , (40) 

d R = U d R d' R , e R = U R e' R , (41) 

where the prime denotes mass eigenstates. Since the up-quark matrix is diagonal, Uf 
coincides with the CKM matrix. 

We can now express the proton decay operators of Eq. (|23|) in term of mass eigen- 
states. 



Aff = 6 °/37 e °k ( U R) kj U %i ( d P,m ( U t) im u l,3 ~ u P,i i U t) fl d' lt l) 



.V 



+ <« (K) lk u%, (Ul) kj e' 3 - d> Jim (Ut) im {Ul) kj */)] + h.c. (42) 

Note that for the orbifold construction, where only the first generation weak eigenstates 
couple to the X bosons, the effective operators read instead 

+ (43) 
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Let us analyze the mixing pattern for the down quarks in the orbifold model. Since 
the mass matrices of down quarks and charged leptons both have the form Eq. (J3J), 
m d ~ m e ~ m, and £/£ have the same structure but, in general, coefficients 0(1) 
will be different. For /Ui,//2 <C ^3, we obtain for the right-handed down-type quarks and 
left-handed charged leptons, 



v 4 v 3 



( M 2 

M12 

Mi 
M12 



V 







M^nzMz-^Mi) 
p, 3 M M12 

M 2 (ji3M3- l j, 3 M i ) 
p 3 M M 12 

P-3 M 

V3 M12 

P-3 M 



Ml( At 3M4+/i3Af 3 ) 

M3 M 2 
M2(m3M 4 +M3M 3 ) 
M3 M 2 

M3 M 3 M 4 -M3 ( M 2 +M 2 +M 2 ) 
P3 M 2 

JL3(M 2 +M 2 +M 2 )-V3M3M 4 

M 2 



Mi \ 

M 

M2 
M 

Ms 
M 



Ma 
M 



(44) 



/ 



up to the two-dimensional mixing matrix for the second and third generation discussed 
in Section 2. 

The rotation matrices of the left-handed down quarks and right-handed leptons are 
obtained by diagonalizing m T m, which leads to small mixing angles, with 



Uf. = V n 



u 



R ■ 



(45) 



The unitary matrices C7£, Ur, and determine the coefficients of the proton decay 
operators in Eq. (f4*3*|) . 

To make a comparison with ordinary 4D GUT models, we consider the flavour struc- 
ture of two SU(5) models described in Refs. [33, 34]. These models make use of the 
Froggatt-Nielsen mechanism [35] where a global U(1) F flavour symmetry is broken spon- 
taneously by the VEV of gauge singlet field $ at a high scale. Then the Yukawa couplings 
arise from the non-renormalizable operators, 

9a f V ) • (46) 



h 



ij 



A 



Here, are couplings O (1) and Qi are the U(1) F charges of the various fermions. 
Particularly interesting is the case with a 'lopsided' family structure, where the chiral 
charges are different for 5* and 10 of the same family. The two examples [33, 34] with 
phenomenologically allowed lopsided charges are given in Tabled 

Note that in these models the large neutrino mixing is explained by a large mixing 
of 5*-plets, which is analogous to the large mixing of lepton doublets and right-handed 
down quarks in the 6D model described above. Contrary to the 6D model, this does 
not determine the U(1) F charges of the right-handed neutrinos. For proper choices these 
models also lead to successful baryogenesis via leptogenesis [36]. 
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Qf 


10 3 


10 2 


10i 


°3 


°2 


51 


model A 





1 


2 


a 


a 


a + 1 


model B 





3 


5 








2 



Table 1: U(1) F charges of the SU(5) fields; a = 0, 1. 



The charge assignments determine the structure of the Yukawa matrices. In model 
A [33], corresponding to the semi-anarchical model of [34], the couplings for down quarks 
and charged leptons read 



/ e 3 e 2 e 2 



\ 



\ 



(47) 



1 / 



where the parameter e = ($)/A ~ 1/17 controls the flavour mixing. Diagonalisation of 
the Yukawa matrices yields 



/l e e 2 \ 



U L = V CKM ~ 



1 e 



1 / 



/ 1 6 6 \ 

e 1 1 
\e 1 I/ 



(48) 



In model B, which corresponds to the hierarchical Hjj model of [34], the structure is 
similar, while the small parameter is instead A ~ 0.35, such that A 2 ~ e. 



4.2 Decay rates and branching ratios 

To calculate the decay rates, we have to evaluate the hadron matrix elements {PS \0\p), 
which describe the transition from the proton via the three-quark operator O to a pseudo 
scalar meson. The various matrix elements are calculated from the basic element 

aP L u p = e a/ 3 7 <^0 

with the aid of chiral perturbation theory [37,38]; u v denotes the proton spinor. The decay 
rates for the different channels are given in Table El Here m p , m n , mx and m v denote 
the masses of proton, pion, kaon and eta, respectively, and f n is the pion decay constant; 
m B = 1.15 GeV is an average baryon mass according to contributions from diagrams 
with virtual E and A; D = 0.80 and F = 0.46 are the symmetric and antisymmetric 
SU(3) reduced matrix elements for the axial- vector current [39]. The matrix element 
a is evaluated by means of lattice QCD simulations; its absolute value varies in the 
range (0.003 — 0.03) GeV 3 . We will choose \a\ = 0.01 GeV 3 (see the recent discussion in 
Ref. [40]). 



u 



p 



(49) 
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T(p -> e+vr ) 

r(p -> UjTT + ) 

r(p - ^+) 

r(p -> e+77) 



327rm3/2 G ^ ^ 

{m2p ~ m ' ±)2 a 2 A 2 GUl + D + FfC 2 



udue-j 



(m 2 -m 2 K0 ) 2 



r 2 

2 

uddv 



m r 



2 

USUCj 



a 2 A 2 G 2 G 



[m 2 p - m 2 K± ) 2 



jm 2 m 2 ) 2 

32xmu; a A 



udsvi 



G 



1 + D-3F 
V6 



C 



udue ; 



Table 2: Partial widths of proton decay channels [41]. 

In the decay rates listed in Table [2] the coupling constant Gq = olIiM^) 2 is given 
by Eqs. (pjDj) , (pTTjl and (JSHJ)- The operators have to be evolved from the GUT scale down 
to the hadronic scale, which is described by the factor A = A SB ■ A LB . It contains both 
a short-distance contribution A SD = 2.37, for the evolution from the GUT scale to the 
SUSY-breaking scale, and a long-distance contribution A LD = 1.43, for the evolution 
from the SUSY-breaking scale to 1 GeV [13]. The effect of lepton masses is neglected. 

The quark and lepton mixing patterns discussed above fix the coefficients C 2 , km in 
the decay rates. As an example, for the process p — > e + 7i°, we obtain from Eq. (l4"3j) . 

4[(^)n(^) 11 ] 2 +[(C/l) 11 (C/2) 11 ] 2 . (50) 



G 2 

udue 



The coefficients for the other processes can be read off analogously. The decay rates in 
Table El have the same form as the decay rates determined by dimension-5 operators. 
The difference lies in the coefficients Ci^ m and the coupling constant Gq- 

We now start with the simplest case of our orbifold model, with U% = U[, and with 
degenerate masses M and /I 3 = // 3 , which we denote I. The mixing matrix for 

right-handed down quarks and left-handed charged leptons, Eq. (J44j), is then simply given 
by 

/ 



U 



R 



1 

V2 





1 

2 




1 





1 

2 


1 
2 





1 

V2 


1 

2 


1 
2 





1 

V2 


1 

2 


U 



(51) 
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decay channel 




Branching 


Ratios [%] 




6D SO(10) 


SU(5) x U(1) F 




case 


I case II 


models A & B 


e+7r° 


75 


71 


54 




4 


5 


<1 


U7l + 


19 


23 


27 


e + K° 


1 


1 


<1 


fi + K° 


<1 


<1 


18 


vK + 


<1 


<1 


<1 


e + r] 


<1 


<1 


<1 




<1 


<1 


<1 



Table 3: Resulting branching ratios and comparision with SU(5) x U(1) F . 



thus the state d\ has no strange-component. For the current-current operators we then 
obtain from Eq. 



eff ~ (Mf) 2 60/87 



2 Kfd eC M « ^M 7 + ^^«^« 7 e + 2 V ud V us /i c u c a dp w 7 
+ 2 KdKs < S/3 w 7 + 2 V^ s Jfv^spu^ 



+ h.c. , (52) 



where the fermions are now mass eigenstates. From this equation we can read off the 
coefficients of the various coefficients Cijkm appearing in the decay rates (cf. Table EJ), 

r 2 = a t/ 4 4- - r 2 = 4 v 4 r 2 =r 2 =av 2 v 2 

^udue ^ ' ud 1 ^ ) ^usufj, ^ v us ) ^uduu ^usue ^ v us v ud 5 

^uddi/ = 2 ' Cnds^ — Kxs , C us du = , (53) 

where we have used Y^j=\ (^zOij (^l)ij = x - 

The numerical results for the branching ratios are listed in Table 01 Note that the 
effect of the derivative operators is negligible for c 5 = c 6 = 1. For the listed branching 
ratios the corrections are less than 3%. 

To compare the branching ratios of the 6D model with those of the two 4D GUT 
models described above, we assume that some mechanism suppresses or avoids the proton 
decay arising from dimension-5 operators. The coefficients C 2 ^ kl can then be derived from 
Eq. (|4*2*|) using the mixing matrices given in Eq. (j48)l . For model A, they read 
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C 2 , ~ C 2 ~ 4 C 2 , ~ C 2 ~ 4 e 2 

^ udue ^usu^i ' ^ udu^i ^usue v ) 

^uddu — 2 ' C u dsu — ~^j= , C us d v — 2 e . (54) 

In model B, e is replaced by A 2 . The differences between the branching ratios of the two 
models are not significant. 

The difference between the 6D SO(10) and the 4D SU(5) models is most noticeable 
in the channel p — > ^i + K°. This is due to the absence of second and third generation 
weak eigenstates in the current-current operators and the vanishing (12)-component in 
Uft and £/£ in the case of the 6D model. Hence, the decay p — > [i + K° is doubly Cabibbo 
suppressed. This effect is a direct consequence of the localization of the 'first generation' 
to the Georgi-Glashow brane. 

Let us now consider the general case, where the M^ d,e ^ are not degenerate, and where 
/i3 and /4 d ' e ^ differ as well. From Eq. (|4^|l we see that the strange component in d\ does 
not vanish anymore, but it is smaller than the bottom component. We have studied 
several cases whose results agree remarkably well. As an illustration, consider the case 
where /xf = 2/x 3 and jl^ = 3fi3, with non-degenerate heavy masses Mf : : M| : Mf = 
I : 73 : 73 : 1 and ^[ : ^2 : M$ : M| = \ : ± : 1 : \ (case II). The branching ratios 
are listed in Table El the differences between the two cases are indeed small. 

The most striking difference is the decay channel p — > fi + K°, which is suppressed by 
about two orders of magnitude in the 6D model with respect to 4D models. It is therefore 
important to determine an upper limit for this channel in the 6D model. Varying the 
mass parameters in the range Mj/M = 0.1 — 1 and Jlf e / fi-s = 0.1 — 10, we find 

Finally, a limit on the compactification scale can be derived from the decay width 
of the dominant channel p — > e + ir°. Neglecting the suppressed contributions from the 
derivative operators, we obtain the analytic expression 
K-m 2 ) 2 fl + D + F^ 



r < p - e v » - -ii^JT aA \rir- 1 (56) 



* 2 v i + 2.3 



16 Mi V V M, 



M d2 Ml 2 

Wud + 



Mf 2 + M$ 2 Mf + Mf 



/9 x 10 15 GeV\ / a \2 , N _ x , 

- 77 d 5.3xl0 33 yrs) . 57 

V M c J Vo.01 GeV 3 / v y ' y ' 

With M* = 10 17 GeV and Mf; 2 e = 0(1), the experimental limit r > 5.3 x 10 33 yields 

M c > M c min ~9x 10 15 GeV, which is very close to the 4D GUT scale. The lower bound 
corresponds to M*/M c min = 12 and (Mf) 2 /(Mf Q ) 2 = 3.72. 
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5 Conclusions 



We have studied proton decay in a 6D SO(10) orbifold GUT model, which is deter- 
mined by dimension-6 operators. The characteristic features of the model are the differ- 
ent breakings of the SO(10) symmetry at different points in the extra dimensions and 
the associated localization of some quarks and leptons. We find that, like in 5D orb- 
ifold GUTs [22], the proton decay rate is enhanced and the branching ratios are strongly 
affected by the quark-lepton 'geography'. 

The summation over Kaluza-Klein towers depends logarithmically on the cutoff scale 
M* in 6D, contrary to 5D, where the sum is finite. Identifying the cutoff with the 6D 
Planck mass, i.e. M* ~ 10 17 GeV, the SuperKamiokande bound on the proton lifetime 
leads to the lower bound on the compactification scale M c > 9 x 10 15 GeV. On the other 
hand, the approximate unification of gauge couplings suggests M c ~ gut — 2 x 10 16 GeV. 
This yields the proton lifetime r(p — * e + 7r°) ~ 1 x 10 35 yrs which, remarkably, lies within 
the reach of the next generation of large volume detectors! 

The peculiar flavour structure of our orbifold GUT model leads to characteristic 
signatures in the branching ratios of proton decay, in particular the strong suppression 
of the mode p — > ji + K° compared to the predictions of 4D models. The reason is the 
higher-dimensional quark-lepton 'geography' and the related non-universal couplings of 
GUT bosons to fermions. Such a pattern can be tested already with a handful of events! 
Hence, the discovery of proton decay may not only confirm the most striking prediction 
of grand unification, but it might also reveal its higher dimensional origin. 
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J. Ellis, and A. Romanino, C. Scrucca and the other participants of the CERN Phen Club 
on 26/2/2004 for useful discussions and suggestions. The work of D.E.C. was supported 
by Fundagao para a Ciencia e a Tecnologia under the grant SFRH/BPD/1598/2000. 
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